Universal adiabatic dynamics in the vicinity of a quantum critical point. 
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We study temporal behavior of a quantum system under a slow external perturbation, which drives 
the system across a second order quantum phase transition. It is shown that despite the conventional 
adiabaticity conditions are always violated near the critical point, the number of created excitations 
still goes to zero in the limit of infinitesimally slow variation of the tuning parameter. It scales with 
the adiabaticity parameter as a power related to the critical exponents z and v characterizing the 
phase transition. We support general arguments by direct calculations for the Boson Hubbard and 
the transverse field Ising models. 



Quantum phase transitions have attracted a lot of the- 
oretical and experimental attention in recent decades, 
see for example Ref. They are driven entirely by 
quantum fluctuations and occur at zero temperature. In 
this paper we will be interested in second order tran- 
sitions, which are characterized by universal properties 
near the critical point. Usually these properties can be 
revealed experimentally by measuring various correlation 
functions. Since the relaxation time in most of conven- 
tional condensed matter systems is relatively short, only 
equilibrium or steady state regimes are experimentally 
relevant. On the other hand recent progress in the re- 
alization of ultra cold atomic gases^ made it possible 
to study experimentally both equilibrium and strongly 
out of equilibrium properties of the interacting quantum 
systems. Thus observation of the superfluid-to-insulator 
transition^ relied on the reversibility of the phase coher- 
ence after the system was slowly driven to the insulating 
state and then back. In the same experiment another 
resonant feature was observed if the Mott insulator is 
a subject to an external linear potential of a particu- 
lar strength. This feature was later interpreted later as 
an Ising-like quantum phase transition between normal 
and dipolar states^. The other big advantage of atomic 
systems is that the parameters governing the transition 
can be tuned continuously during a single experiment, 
so that, for example, it is possible to cross a quantum 
critical point in a real time. 

Let us consider now a specific situation, where some 
system was initially in the ground state. Then a tun- 
ing parameter was slowly changed to drive it through a 
critical point. From general principles we know that the 
system should remain in the ground state as long as it 
is protected by the gap from the excitations. On the 
other hand the gap vanishes right at the critical point 
so the adiabaticity conditions can never be satisfied in 
the vicinity of the phase transition. The slower the pa- 
rameter changes the more time the system spends near 
the critical point, but on the other hand the less the 
interval where the adiabaticity is violated. The compe- 
tition between these two processes determines the total 
amount of excitations in the system. Here we show that 
the number of excited states decreases as a power law of 
the tuning rate. Because of the universality and scaling, 



below a certain dimension, which we identify later, this 
power is determined by the critical exponents z and v 
characterizing the transition. So measuring the number 
of excitations as a function of the tuning rate one can ob- 
tain the information about the critical properties of the 
phase transition. We give a general argument for the par- 
ticular scaling form and consider two specific examples of 
phase transitions occurring within Boson Hubbard- and 
transverse field Ising models, which confirm this scaling. 

Let us start from a general formalism. We assume 
that the system is described by some Hamiltonian H(\), 
which depends on the external parameter A. Without loss 
of generality A = corresponds to the phase boundary, 
so that A > and A < describe different phases of the 
system. Let the set of (many-body) functions <f> r {X) rep- 
resent the eigen-basis of the Hamiltonian H. The wave- 
function of the system can be always expanded in this 
basis: 

V; = 5>p(<)0 P (A). (1) 

p 

We assume that A slowly changes in time: A(i) = 6t, 
where S is the adiabaticity parameter and we took linear 
dependence on time for the sake of convenience. Then 
substituting (JIJ into Schrodinger equation we find: 

i^ +i 5j2a q (t){p\-^\q)=LJ p (X)a p (t), (2) 

where w p (A) is the eigen frequency of the Hamiltonian 
Ti{\). It is convenient to perform a unitary transforma- 
tion: 

a p (t) =h p {t)e- l $ t ^ {x{t))dt = h p {\)e-i^ ul ^ x)dX . (3) 
Combining @ and we derive: 

^ = -E^)(pilk)^ /AK(V) ^ (V))dA '- (4) 

q 

If before the evolution the system was in the ground state 
|0) then a single term dominates the sum in (0J. The 
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relative number of the excited states is thus given by: 



E 

p 



dX(p\ — \0)e 



■ J x (u lp (\')-uj (\'))dX' 



(5) 

where the prime over the sum implies that the summation 
is taken only over the excited states. It is important to 
emphasize that across the second order phase transition, 
which we consider in this paper, the basis wave functions 
change continuously with A. 

Let us assume that we deal with a uniform d- 
dimensional system. This assumption is not necessary, 
but it is the case for the most known systems undergoing 
a quantum phase transition. We also assume that there 
is a single branch of excitations characterized by the gap 
A and some dispersion. Since both A and A become 
zero at the phase boundary then, near the critical point, 
we must have A cx |A| Z1/ [lj with z and v being critical 
exponents. In the momentum space (0 reduces to: 



d d k 
(2n) d 



d\(k\ — |o)e^ A ^ A ')-^ A '» dV 



From general scaling arguments we can write: 
uj k -cj = AF(A/F) = V F(\™ /k z ), 



(0) 



(7) 



where F (F) is some undefined function satisfying 
F(x) cx 1/x for large x, and z is the dynamic critical 
exponent. Similarly we can argue that: 



-G(A/fc 2 ) — 10) = —j^—G{\ zv /k z ), 

(8) 

where G (G) is another scaling function satisfying G(0) = 
const. Having these scaling forms in mind we can do the 
following substitutions in (0: A = k 1 ^^, k = S zij + 1 r]. It 
is easy to see that if the momentum integral in (0 can 
be extended to infinity then: 



(9) 



where C is a nonuniversal constant which depends on 
the details of the transition. The condition allowing to 
send the upper cutoff to infinity is d < d c = 1z(zv + 1), 
where d c is the upper critical dimension for this prob- 
lem. Note that v and z can depend on d themselves. For 
d > d c , the main contribution to the excitations comes 
from the high momentum states. In this would 
still vanish at S — > 0, but the universality will disap- 
pear as excitations with high momenta (k z 3> A) will 
dominate the transitions. The result (0 is quite remark- 
able, it shows that if 5 — ► and v is a finite number 
greater then zero, the transitions to the excited states 
are suppressed and the adiabatic limit still holds. Wc 
want to emphasize, the adiabaticity is understood in a 
sense that the density, not the total number, of excita- 
tions is much smaller than one. Strictly speaking in the 



true adiabatic limit there are should be no excitations 
and the system must remain in the ground state. How- 
ever to achieve this, it is necessary to scale 5 as inverse 
power of the system size^, which is virtually impossible 
to do in large systems. There are two limits v — > and 
v — > oo in 0) which require special attention. The first 
one is trivial since it corresponds to the absence of the 
phase transition since the gap always remains finite ex- 
cept for a very narrow interval around A = 0. Indeed, a 
more careful analysis shows that the constant C in is 
proportional to v 2 . The opposite limit v — * oo is more in- 
teresting since it corresponds to e.g. Kosterlits-Thouless 
(KT) transition 6 , which has many realizations in 1 + 1 di- 
mensional quantum systems. Thus if the precise scaling 
form is A cx c~ ^" then 



5* hf 



■(s- 1 ). 



(10) 



This expression acquires extra logarithmic corrections as 
compared to 0. In particular, for the KT transition 
r = 1/2 and z = 1 so that (|10fl reduces to 



KT 



cx^ln 3 ^- 2 ^- 1 ). 



(11) 



There are no logarithmic corrections in two dimensions. 
However, the only physically relevant case where the KT 
transition can occur in a quantum system at zero tem- 
perature corresponds to d = 1 . 

Qualitatively one can interpret |5| in a simple way. 
The transitions to the excited states occur when the adi- 
abaticity conditions break down, i.e. when dl ^ t A > A. 
From this one immediately finds that the time interval 
when the transitions take place scales as: t ~ 8~ . 
The typical gap at this time scale is 



A - {8t) zv - Stt, 



(12) 



which amounts to the available phase space f2 ~ k d ~ 
A~ ~ d^TT. Now if we use the anzats that dA(k\d/dA\0) 
is a scale independent quantity (see ©), then we im- 
mediately come to the conclusion that this phase space 
determines the number of excited states n ex so that we 
come to (0. This simple derivation above, in fact, docs 
not rely on the spatial homogeneity of the system. The 
only information we need is the density of states of exci- 
tations p(e) at the energy scale determined by l|12(l . So 
in a general case instead of 10 we get: 



oc S^ 1 p(S^+ 1 ). 



(13) 



Notice that 10 contains only two critical exponents v 
and z. So measuring the dependence n ex (S) and knowing 
one of the exponents, say z, one can immediately deduce 
other. 

Let us apply these ideas to the superfluid-to-insulator 
transition in a system of interacting bosons in a d- 
dimensional lattice at commensurate fillingiit To de- 
scribe the excitations near the critical point we adopt a 
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mean-field Hamiltonian derived by Altman and Auerbach 
in Refs. 00: 

H = 2dJN^2< (2 u cos 9 — cos 29 + ek cos 2 9) b[ k &i,k 



l-£k 



cos 2 9(b\ ik ^_ k + 6i, k 6i,-k) 



(2u cos 2 - + sin 2 9 - cos 2 - + e k cos 2 -) &J, k 6 



— — COS 2 -(4,k 6 2 ,-k + &2,k&2,-k) 



2> "2,U° 2 < k 

(14) 



Here J is the tunneling constant, N is the mean number 
of bosons per site, which we assume to be a large integer 
for the sake of simplicity, e k = wg Yls ^ ~ e ^ S > ® ls the 
mean field angle characterizing the phase. In particular 
9 = corresponds to the Mott phase, while in the su- 
perfluid regime cos 9 w u. The dimensionless interaction 
u = U/(4:JdN) is defined so that the transition occurs at 
u = 1. 

In the Mott side of the transition, u > 1, both branches 
are degenerate and the Hamiltonian (|14fl can be readily 
diagonalized via the Bogoliubov's transformation: 

/3m,k = cosh </> mi k fo m ,k - sinh <p m ^ b] n _ k , (15) 
where it becomes 



W = I>m,k/4 ik /3m,k. (16) 
m.k 

The eigenfrequencies w m . k and the angle 4> m ,^ read: 



-V-.k 4,1-iXaJ- — - + — « AdJNy/X + e k , (17) 



u u 



tanh 2</>i,2. k = ± 



l-e k 



2A + l + e k ' 



(18) 



We have chosen A = u — 1 to be the parameter governing 
the phase transition. Given Hamiltonian (|16fl and trans- 
formations i|15|l it is easy to write down the ground state 
wavefunction: 



|0) = Y\ cosh(/> mik e 

m,k 



tanh (bm ib 



k»L.k»L,-h|Vac), (19) 



where | Vac) is the state with no b particles. It is a simple 
exercise to check that (p\d/dX\0) is nonzero only when 
two particles with opposite momenta are excited, i.e. 



|p) = |m,k,-k) = 6t k fit | 0) . 



(20) 



Then it can be verified that 
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(m,k, k l aA lO)-T 2(1 + A)(A k +£k) - . 2 . A 

(21) 



where we used the approximation that both A and e k are 
small near the phase transition. Note that H21|) satisfies 
the general scaling (jHJ with the exponent v = 1/2, in the 
same way the dispersion relation H17|) agrees with J7|). 

A similar analysis can be performed on the superfluid 
side. The Hamiltonian l|14|) gives now two branches cor- 
responding to the amplitude and the phase modes: 



wi,k ~ 4dJNyJ-\ + e k , w 2 , k w AdJNy/^, (22) 

which are characterized by the following angles of the 
transformation l)15[l: 

♦ ^ 1 - £ k + 2A 2 - 2e k + A 

tanh20i !k « — — , tanh20 2 .k « tt— : — r- 

1 + Ek ^ + *£k + A 

Note that the parameter A is negative on the superfluid 
side. The matrix elements for these two modes are: 



25b)' < 2 ' k -- k i|;i ^ 1 
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The excitations of the phase modes are suppressed as 
compared to the amplitude ones. This can be also ex- 
pected on the physical grounds, i.e. by changing the 
parameter A or equivalently u we change the mass or 
the gap of the amplitude mode thus exciting it, how- 
ever there is no such a coupling mechanism for the phase 
mode. Clearly the total number of the excited phase os- 
cillations is not determined by the properties of the criti- 
cal point and thus is not universal. This number scales as 
n ex tx S d , i.e. vanishes much faster with S than the num- 
ber of excitations to the gapped mode. Besides, a typical 
experiment would use the phase contrast as a measure 
of superfluidity 3 ^, which is not strongly affected by low 
momentum phase excitations. Keeping this in mind we 
calculate explicitly only the number of particle pairs lost 
to the mode 1, which is gapped on both sides of the tran- 
sition. Experimentally this number can be detected by 
getting first from the superfluid to the Mott insulator and 
then returning back to the superfluid regime and measur- 
ing the loss of the phase contrast, or by measuring the 
number of created particle-hole pairs in the insulating 
state. Performing the integration in (JSJ we find that in 
one, two and three dimensions the number of excitations 
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a 0.348 ( 
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3D _ 

ex 


a 0.010 ( 


f S 
v JN 



2D 



0.059 
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In 



(25) 



respectively. We note that since v = 1/2 within this 
meanfield treatment and z = 1 the upper critical dimen- 
sion is 2z(zv + 1) = 3 so the scaling J§J) is valid in one 
and two dimensions. However, this model has an addi- 
tional symmetry, giving the same prefactors of the gap 
and wavefunction dependence on A in both superfluid and 
insulating phases (compare l|17|) with l|22|l and (|21|) with 
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(12411 1. This symmetry shifts the upper critical dimen- 
sion to d = 4, so the results remain universal in all three 
spatial dimensions. As anticipated, expressions (|25[) are 
consistent with J??}. 

The correct description of the superfluid-to-insulator 
transition gives exponents different from the mean-field 
values used above. The commensurate transition lies in 
the universality class of the XY model in the d + 1 di- 
mensions with z = 1 and v — 0.5 for d = 3, v « 0.67 for 
d = 2 and v — oo in d = 1— (more precisely the uni- 
versality class in the latter case is of the KT transition 
with A oc exp(— b/\/~X)). So that for these special points 
equation 10 reduces to: 



7Z„™ OC i-> 



,2D 



oc 5' 



0.80 



n^ocJ. (26) 



In a generic point of the superfluid-insulator transition 
corresponding to the non-commensurate filling zv = l 7 
and since z > 1 the upper critical dimension d c is always 
larger than 3. So © reduces to: n ex oc 6* u . Thus mea- 
suring n ex (S) it is possible to observe the exponent v. 
Unfortunately except for 3D, where v — 1/2 and hence 
n ex oc <5 3 / 4 , the precise value of v (and z) is not fixed 
but rather depends on the point where the transition oc- 

1 9 

curs**. 

Another example we consider here is the transverse 
field Ising model 1 , which is described by the Hamiltonian 



(27) 



where a x and er 2 are the Pauli matrices. The dimen- 
sionless coupling constant g drives the system through 
a critical point, which occurs at g c = 1 and which 
is characterized by the critical exponents z = v = 1. 
Using the Jordan- Wigncr transformation one can show 
that H27|) maps to the model of free spinless fermions 
with the Hamiltonian 

Hi = - ^ cjc i+ i + ct +1 Cj + ct c ] +1 + c j+1 Cj - 2gc\c j , 

(28) 

which in turn can be diagonalized by the Bogoliubov's 
transformation: 



c k = cos(0 fe /2)7 fc +isin(0 fc /2)7! 



(29) 



Here Ck is the Fourier transform of Cj and the angle 9k is 
given byi 



tan 6]~ 



sin(fe) 
cos(fc) — 



(30) 



In the diagonal form the Hamiltonian l|28|l reads 
Hi = ^e k ili kl 



(31) 



where eu = 2yl + g 2 — 2gcosk. The ground state wave- 
function, which is the vacuum of l|31|l reads 

|fi) = Y[(cos(9 k /2) + ism(9 k /2)clcl k )\0) (32) 



where |0) is the state with no c-fermions. The excited 
states have the form of 7L7L • ■ The natural 

choice of the tuning parameter A is A = g — 1 which is 
proportional to the energy gap A. It is straightforward 
to verify that 



0_ 



sin k 



pl-llM, (33) 



g 2 — 2g cos 



which again corresponds to two particle excitations. This 
expression is consistent with © with G(0) = |. Now 
using we immediately find: 



(34) 



which agrees with the general formula © given that d = 
v = z = 1. 

In conclusion, we showed that if the system, origi- 
nally in the ground state, is slowly driven through a 
quantum critical point, the number of excited states per 
unit volume goes to zero as a power law of the tuning 
rate. The exponent is universal and is determined by 
the critical properties of the transition if the dimension 
is smaller then d cr — 2z(zv+ 1). We provided some gen- 
eral arguments and performed explicit calculations for 
the superfluid-to-insulator transition within the Boson 
Hubbard model and for the quantum phase transition in 
the transverse field Ising model. 

Recently, two other papers appeared, which addressed 
a similar issue of the number of created defects for a spe- 
cific case of a transverse field Ising modeii&it. In partic- 
ular, the authors got the same scaling as in Eq. (34) but 
with a slightly smaller numerical prefactor. The discrep- 
ancy comes from a more accurate treatment of transition 
probabilities within the Landau-Zecner formalism 15 . 
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Materials Research Laboratory via grant DMR-0213805. 
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